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Abstract 

We investigate the electrostatic contribution to the persistence length of 
a semiflexible polymer chain whose segments interact via a screened Debye- 
Hiickel dipolar interaction potential. We derive the expressions for the renor- 
malized persistence length on the level of a l/D-expansion method already 
successfully used in other contexts of polyelectrolye physics. We investigate 
different limiting forms of the renormalized persistence length of the dipo- 
lar chain and show that in general it depends less strongly on the screening 
length than in the context of a monopolar chain. We show that for a dipolar 
chain the electrostatic persistence length in the same regime of the param- 
eter phase space as the original Odijk-Skolnick-Fixman (OSF ) form for 
a monopolar chain depends logarithmically on the screening length rather 
than quadratically. This can be understood solely on the basis of a swifter 
decay of the dipolar interactions with separation compared to the monopolar 
electrostatic interactions. We comment also on the general contribution of 
higher multipoles to the electrostatic renormalization of the bending rigidity. 
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1 Introduction 



Semiflexible polyelectrolytes are ubiquitous in biological context ranging from 
biopolymers such as DNA, filamentous F actin, microtubules, and then all the 
way to molecular aggregates such as micelles or even whole organisms as in bacte- 
rial fd viruses or the tobacco mosaic virus. In all these cases we are dealing with 
objects that at an appropriate scale behave as Euler-Kirchhofhan elastic filaments. 
The most important mechanical characteristics in these molecular systems is the 
persistence length stemming from the bending rigidity of the polymer that can 
exhibit enormous variation in magnitude bracketted by nm on the lower end and 
mm on the upper end. 

The bending rigidity and thus the persistence length is a consequence of short 
range atomic and molecular interactions and is itself a mesoscopic property pp. 
For charged semiflexible polymers the long range nature of the electrostatic in- 
teractions modifies the value of the persistence length as is well known from the 
seminal work of Odijk, Skolnick and Fixman (OSF) (2j. Recent detailed critical 
assessment of the OSF conjecture jSHHEj confirmed the universality of the depen- 
dence of the persistence length on the parameters of the electrostatic interaction. 
It appears that the OSF behavior characterized by the inverse dependence of the 
persistence length on the ionic strength of the bathing medium is quite robust. 
One of the main ingredients of the OSF ideology is a complete lack of any atomic 
or molecular specificity in the interaction between charge segments along the poly- 
electrolyte. By assumption the mobile charges in solution and the fixed charges on 
the polyelectrolyte backbone interact only via generic (screened) Coulomb inter- 
action. Though this assumption seems to be reasonable it would be appropriate to 
investigate what kind of behavior of persistence length is conferred by the effects of 
strong specificity as in e.g. the case of specific ion adsorption. Recent experiments 
also suggest that reevaluation of the basic assumptions in our understanding of 
counterion - backbone interactions in polyelectrolytes would be highly desirable 

In specific ion adsorption mobile charges from the bathing solution and fixed 
charges along the polyelectrolyte can combine, leading to the emergence of higher 
multipoles along the polyelectrolyte chain, the first one being a dipole stemming 
from the association of a negative fixed charge and a specifically adsorbed mobile 
positive charge from the bathing solution. The same situation though for different 
reasons, specifically due to the low dielectric constant envoronment, could be ob- 
tained also in the context of ionomers [Zj. Similar consideration could apply also 
to electrorheological or ferrofluids in external fields. 

Muthukumar [Hj was the first to realize the importance and the extent of mod- 
ifications wrought by the first higher multipole i.e. the dipole, on the behavior of 
a flexible polyelectrolyte described on the level of the Edwards Hamiltonian. He 
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discovered the formation of localized aggregated structures along the chain that 
dominate the statistical behavior of the flexible polyelectrolyte chain. They are 
characterized by a different scaling of the size of the chain with respect to its length 
and depend continuously on the parameters of the dipolar interaction potential. 

Chain flexibility as formalized by the Edwards Hamiltonian [5] is essential for 
emergence of this type of aggregated structures. Now assume a dipolar semiflex- 
ible chain, described as an Euler-Kirchhoffian filament with dipolar intersegment 
interactions. On top of that assume that the length of the chain is the largest 
length in the problem. This constitutes a complementary limit to be contrasted 
with the Muthukumar calculation. The analysis provided by Muthukumar and the 
one detailed below should thus brackett the behavior of any real polyelectrolyte 
chain with dipolar charges along its length. 

The problem with dipolar interactions is that they depend on the local ori- 
entation of the polyelectrolyte segments. This complicates the evaluation of the 
partition function of the semiflexible chain in an essentail manner. In order to make 
this complicated problem tractable, we will invoke the 1/D-expansion method jIT] 
that has already been successfully applied to the problem of persistence length 
renormalization of a semiflexible polyelectrolyte chain ^2] with monopolar charges. 
This will allow us to explicitly evaluate the electrostatic contribution to the per- 
sistence length of a dipolar chain. As a point of departure we will assume that 
the dipolar charges along the chain interact via a screened Coulombic interaction 
potential, the range of which is characterized by the Debye screening length [T3j . 
We will derive the complete dependence of the electrostatic persistence length on 
the parameters of the dipolar interaction potential and show that it leads to a 
different behavior regarding its dependence on the ionic strength of the bathing 
solution then the standard OSF result. We find basically two different regimes for 
the bahavior of the electrostatic persistence length of a dipolar chain 

• for small but non- vanishing values of the screening length Ad the electrostatic 
contribution to the persistence length behaves as 

Lp = Lp + const, log Ad 

where Lp^ 1 is the bare value of the persistence length. This expression would 
superseed the OSF X 2 D result valid for monopolar interactions. Obviously 
the dipolar electrostatic renormalization of the persistence length depends 
more gently on the screening length. 

• for large values of the screening length Ad the electrostatic contribution to 
the persistence length behaves as 

LpP = const. A o 13 , 
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where /3 equals either 3/4 or 3, depending on the strength of the dipolar in- 
teraction. This result should be compared with various variational estimates 
of a sub-OSF behavior for a monopolar chain that lead to (3 in the vicinity 
of 1. 

These two regimes stemm from two different limiting forms of a single equation 
giving the renormallized value of the persistence length as a function of the param- 
eters describing the dipolar interactins along the chain. We also present complete 
numerical solutions of this equation in various regions of the parameter space. 

The outline of the paper is as follows. In Section 1 we first rewrite the Hamil- 
tonian of a semiflexible polymer chain with screened dipolar interactions, assumed 
to be composed of the Euler-Kirchhoff elastic energy and the interaction energy, in 
a form that allows for a straightforward application of the l/D-expansion ansatz. 
We assume that the dipoles of the polymer segments are oriented along the local 
tangent vectors of the chain. This effective Hamiltonian captures the elasticity of 
the chain, the inextensibility of the chain and the fact that electrostatic interac- 
tions depend on the position of the interacting segments as well as their orientation. 
Section 2 introduces all the important approximations in order to make the evalu- 
ation of the partition function of the chain tractable. A diagonalization ansatz is 
introduced for the orientational part of the Hamiltonian, a global inextensibility 
constraint is substituted for the local one and a saddle point evaluation is intro- 
duced for all the auxiliary fields that do not enter the Hamiltonian on a quadratic 
level. In section 3 the saddle point equations are solved explicitly by expanding 
all the Fourier components of different auxiliary fields to the fourth order. This is 
consistent with the semiflexible ansatz for the configuraional part of the Hamilto- 
nian of the chain. Section 4 introduces explicit equations for the bending rigidity 
renormalization that follow from the saddle point equations of the previous sec- 
tion. These equations are solved numerically in the last section. Different limiting 
forms of the numeric solutions for limited regions of the parameter phase space are 
derived also analytically. We take a critical look at the results derived in this work 
and comment on the implied limitations of their validity and their connection with 
previous work on the electrostatic renormalization of the rigidity of semiflexible 
polymers. 

2 The effective Hamiltonian 

We investigate a semiflexible polymer chain with dipolar charges along the chain. 
The dipoles can be either structural or they can be a consequence of the specific 
adsorption of e.g. positive mobile charges from the bathing solution onto the fixed 
negative charges on the polyelectrolyte backbone. The total interaction energy of 
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the chain is given by 



\ I [ V(r(s),r(s'))dsds' 
1 Jo Jo 



(1) 



where 

V(r(s),r(s')) = 4 J 7 r '(t) - r(loi 5 ^ " r(g ' )|2(1 + ^ " r{s ' )l)p{s) ' p(gV 

-(3 + 3«|r(s) - r(s')| + («|r(s) - r( S ')|) 2 )(p(s) • (r(s) - r(s')))(p(s') • (r(s) - r(s')) 

with |r(s) — r(s')| the separation between two segments located at r(s) and r(s') 
along the chain, s is the arclength of the chain, and p(s) and p(s') are dipoles per 
unit length located at s and s'. The above form of the screened dipolar interaction 
follows straightforwardly from the second order multipole expansion of the screened 
Coulomb kernel [H] and reduces to the usual form of the dipolar interaction in 
the lim it of no screening. On the Debye - Hiickel level k = A^ 1 where Xd is the 
Debye screening length. We furthermore assume that the dipole per unit length 
along the chain is given by 

p(s) = p d s r(s) (3) 

where po is the strength of the dipolar moment per unit length of the segment 
and d s r(s) = r(s) is the unit tangent vector. It is thus assumed that the dipoles 
point along the chain, and that the component of the dipoles perpendicular to the 
local axis of the chain is averaged out to zero. Or model is thus "perpendicular" 
to ther case investigated in the context of protein folding jT^j. The interactions 
along the chain are therefore described with a generic form of segment - segment 
interactions that can be cast into the form 

V(r(s), r (s')) = V R (r(s) - r (s'))d s r(s)dMs') - 

-V P (r(s) - r(s')) (d s r(s) ■ (r(s) - r(s'))) (^r(s') ■ (v(s) - r(s'))) , 

(4) 

with 

e -K|r(s)-r(s')| 

V R (r(s) - r(s')) = ^o^^ — -^(l + /t|r(s) -r(s')l), 

V P (v(s) - r(s')) = v i^'t^OP ( 3 + 3 ^ |r(g) " r(g/)l + * 2|r(s) " r ( g/ )| 2 )(^) 

Obiosuly the interaction Eq. 0] depends on the positions of the interacting seg- 
ments as well as on their orientation. This is the fundamental difference between 
monopolar and dipolar interactions. We have defined 

Po _ 7„ m „ fPo 



vo = -r^- = k B T £ B ^ , (6) 
4vree V e o 
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where Is is the Bjerrum length [T3j. The units of vq are energy times length and 
are thus the same as the units of bending modulus of a semiflexible polymer chain. 

The generic form of the interaction Eq. Altogether with the Euler-Kirchhoffian 
elastic conformational part of the mesoscopic free energy jTHj can now be used to 
investigate the statistical behavior of a semiflexible dipolar chain. The mesoscopic 
Hamiltonian of the chain can be written canonically as 



n[r(s)\ = \kcJ\s {d*r(s)) 2 + ± ££ dsds'V(r(s), r (*'))■ 



(7) 



where the local curvature is 



dr( 



d 2 r( 



ds 2 



The bending rigidity Kc is connected with the bare persistence length via Kc = 
ksT Lp' 1 . Since in what follows the chain is assumed to be inextensible an addi- 
tional constraint of d s r(s)d s r(s) = 1 should be taken into account. The partition 
function thus assumes the form 

S = J V[r(s)]U s S 3 ((d s r(s)) 2 - l) exp -(3H\v{s)}. (9) 

There are two difficult problems connected with this partition function: first of 
all we have the constraint of inextensibility that has to be enforced locally for 
every conformation of the chain, and on top of this there is the non-locality of 
the segment-segment interaction potential that couples different segments along 
the chain and depends also on their orientation. We will address both problems 
systematically by applying the method of Lagrange multipliers |T2] . 

To this effect we will introduce two additional terms into the interaction Hamil- 
tonian. The inextensibility constraint is dealt with via the Lagrange multiplier 
A(s), leading to an additional term in the Hamiltonian of the form 

together with an additional trace over A(s) in the expression for the partition 
function. In what follows we will assume that the inextensiblity constraint can 
be implemented globally instead of locally J7j. This automatically means that 
A(s) = A. The non-locality of the interaction potential is more difficult to deal 
with then in the case of a chain without orientationally dependent interactions. 
In order to pave the way for an approximate evaluation of the partition function 
we introduce two new Lagrange multipliers in the form of two tensorial auxiliary 
fields. First of all we define 

B lk (s, s') = (Ti(s) - v t (s')) (r k (s) - r k (s')) (11) 
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and then 

T ik (s, s') = d s ri(s)d s r k (s). (12) 
We notice immediately these two fundamental identities 

Tr B ik (s, s') = (r(s) - r(s')) 2 and Tr T ik (s, s') = d s r(s) ■ <9 s r ( s). (13) 

We thus conclude that Tr Bi k (s,s') is equivalent to B(s,s') introduced in the 
context of orientationaly independent intrachain interactions ^2]- Introducing 
furthermore the following two additional coupling terms 

5H 2 = ^ J J dsds'g ik (s, s') ((i-j(s) - r^s')) (r fc (s) - r k (s')) - B ik (s, s')) 

5H 3 = ~JJdsds'p ik (s,s')(d s ri(s)d s r k (s)-T ik (s,s')). (14) 

we can write the partition function in the following form 

P[r(s)]P[A( S )]P[^ fe ( S )]P[^ fc ( S )]X>[p ife (s)]P[^ fe ( S )] exp -0H[t( 8 )], (15) 



where the effective Hamiltonian 7-^[r(s)] is given by 

H[r(s)\ = H[r(s)\ + 5Hi[r(a)] + 5H 2 [r(s)\ + 5H 3 [v(s)\. (16) 

The new auxiliary fields gi k (s,s') and Pi k (s,s') thus act as tensorial Lagrange 
multipliers setting the constraints (r;(s) — rj(s')) (r k (s) — r k (s')) — Bi k (s,s') = 
and <9 s rj(s)<9 s r fc (s) — T ik (s, s') = 0. In the new variable space the segment-segment 
interaction potential can be obtained in a more transparent form as 

V(r(s),r(s')) = V R (TrB lk (s,s'))TrT tk (s,s')- 

-Vp(Tr B ik (s, s'))B ik (s, s')T ik (s, s') = V (B ik (s, s'), T ik (s, s'( 



In order to arrive at a more transparetn form of the partition function we further- 
more introduce the following two new variables 

\-fc( s ' s ') = >^ S ik 5(s - s') +p ik (s,s') 

gS(s,s') = h(s-s')Jds"(g ik (s,s")+g ik (s»,s'))-g ik (s,s'). (18) 
With these new definitions the effective Hamiltonian can be finally reduced to this 
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fairly complicated form 




J j dsds'2g\ c k '(s,s')ri(s)r k ( 
dsds'p ik (s, s')T ik (s, s') + 



The above form of the Hamiltonian together with Eq. represents at this stage 
an exact expression for the partition function. This is the starting point of different 
approximations introduced below. 

Though the above form of the effective Hamiltonian looks prohibitively com- 
plicated it can in fact be reduced to analytical quadratures provided one devises a 
powerfull enough approximation scheme. It was shown in recent work fF2] ITH| ITH] 
that the 1/D-expansion method, where D here and below is the dimensionality of 
the embedding space, can be fruitfully applied to polymer problems of the above 
type and we will use our experience gained in the context of monopolar interac- 
tions ^2] to tackle also the more complicated case of multipolar interactions as 
implied by the Hamiltonian Eq. E3 



Our rationale for writing the effective Hamiltonian in the form Eq. EI]is that it will 
be shown to be amenable to straightforward approximations leading to a closed 
form evaluation of the partition function. Similar methods have been already used 
successfully in the context of monopolar interactions. First of all we will introduce 
a diagonalization ansatz of the form 



Though this ansatz is not absolutely necessary in order to make progress in the 
evaluation of the partition function Eq. it certainly makes the whole calculation 
manageable and controllable, reducing its overall complexity to a bare minimum. 
The above ansatz clearly implies also 



3 The ansatz 




g(s, s')5 ik 
p(s,s')5 ik . 



(20) 




(21) 
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In order to be consistent we must also introduce the same type of diagonaliza- 
tion ansatz for the fields Bik(s, s') and %k(s, s'), which are conjugate to the fields 
9ik(s, s') and Pik(s, s'). We formulate this as 

B ik {s,s') = jjB(s,s')5 ik 

% k (s,s') = ±T(s,s')5 lk . (22) 

In this way we can write Tr Bi k (s, s') = B(s, s') and Tr T ik (s, s') = T(s, s'). Since 
we are not concerned here with explicit analysis of finite size effects we can assume 
overall that the (s, s') dependencies actually reduce to (s — s'). Similarly to the 
case of orientationally independent intrachain interactions we see that the partition 
function is now quadratic in r(s) and its first and second derivatives. We can thus 
trace over these harmonic degrees of freedom [12J. Expanding the whole effective 
Hamiltonian around a reference configuration r (s) and introducing the Fourier 
components of different fields in the standard way we obtain after tracing out the 
harmonic degrees of freedom the following form of the effective Hamiltonian 

n' = n[ro(s)} + ^f-E Q ^gG(Q)-lxjds + 

~\ J J dsds '9(^ s ') B (^ s ') ~ \ J J dsds'p{s,s')T(s, s') + 

+ dsds'V(B(s,s'),T(s,s')) } (23) 

where the self energy function G{Q) is given by 

Q{Q) = {K C Q A + X^(Q)Q 2 + 2g^(Q)) . (24) 

Here X^(Q) and g <yC \Q) are of course the Fourier components of the fields g <yC \s — 
s') and A^(s — s'). The Hamiltonian of the reference configuration Tq(s) can be 
derived as 

H[r (s)} = l -K c j ds{d 2 Ms)) 2 + ^JJ dsds'X^(s - s') (d s r (s) ■ d s r (s)) + 
+ IJJ dsds'2g^\s - s') (r ( S ) ■ r (s)) . (25) 

All along we disregard the finite (chain) size effects and assume that the chain is 
homogeneous and s, s' dependence is actually equivalent to \s — s'\ dependence. 
This assumption furthermore implies that the length of the chain is the largest 
length in the problem. Thus with this line of reasoning 

9 {C \Q) = jducosQugV(u)=g(Q)-g(Q = 0) 

A^(Q) = [du cos Qu A (c) (u). (26) 
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From here we can also straightforwardly conclude that we have the following iden- 
tities for the functional derivatives 



— = l-cosQ(s-s) 

dg(s,s') 



5\(s,s') 



cosQ(s-s'). (27) 



The next logical step in order to evaluate the partition function Eq. would 
be to trace out the auxiliary fields A(s), gik(s),Bik( s ),Pik( s ),%k( s )- Unfortunately 
these fields enter the effective Hamiltonian non-linearly and they can not be simply 
traced over. Thus we have to introduce a new approximation at this point: instead 
of integrating over the auxiliary fields we will simply evaluate the saddle - point of 
H. with respect to these fields. This is easier and most importantly it is feasible to 
do. It is thus at this juncture that a critical step in the derivation of the partition 
function has to be made that constitutes the essence of the 1 / Z)-expansion method 
as applied to the polymer problem. The details and ramifications of this step have 
been addressed in the previous work 



4 The saddle point equations 

The saddle point of the effective Hamiltonian Eq. 123 w.r.t. the fields X(s,s'), 
B(s, s'), g(s, s'), T(s, s') and p(s, s') can be obtained straightforwardly in a stan- 
dard fashion. For the saddle point in the X(s, s') variable we obtain the following 
relation 

1 = (d s ro(s) • <9 s r (s)) + Dk B T ^^ 1 (g)g 2 , (28) 

Q 

The next saddle point that we evaluate is B(s, s'). It leads to the following equation 

9 W(Q) = /^l-cosQ s) ™ ( ;^» . (29) 

Then follows the saddle point in g(s, s'), giving rise to 

B(s, s') = (r„(s) - r (s')) 2 + 2Dk B T ^ ^ 1 (Q)(1 - cosQ(s - s')). (30) 

Q 

Furthermore the saddle point in the variable T(s, s') leads to 

P{Q) = J dscosQs ^nm>i, (31) 
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and finally the saddle point in p(s, s') gives rise to 

T(s,s') = (d s v (s)-d s r (s)) + Dk B T ^G'^Q 2 cosQ(s - s'). ( 32 ) 

Q 

We now assume that we can retain only terms up to the fourth order in Q for the 
functions \( C '(Q) and g^ c \Q)- This is consistent with the fact that the original non- 
interaction part of the Hamiltonian contains only Euler-Kirchhoffian elasticity and 
is thus at most of fourth order in the Q space. In this way we obtain the following 
expansion 

A««) = A + p(Q) = A + / ^™p> - % / 

(33) 

as well as 

Thus we have basically derived an expansion also for the self-energy Eq. |2U which 
can now be cast into the form 

G(Q) = (A + 8X) Q 2 + (K c + 5K C ) Q 4 + • • ■ = \ {R) Q 2 + K^Q* + .... (35) 

The expansion in Q to the fourth order thus allows us to introduce the renormalized 
values of the parameters A and Kc via 

\( R ) = \ + 5X K [ C R) = K C + 5K C , (36) 

where we introduced the following two abbreviations 



&T(s) J dB(s) 

SKo = -i/<b^2^-A/*^223^M. (37) 

At this point we are in a position to evaluate all the Q integrals in the saddle point 
equations. Thus instead of Eqs. 1281 1501 and 1821 we remain with 

1 = (d s r (s) ■ <9 s r (s)) + Dhl ' 1 



B(s,s) = (r (s) - r Q (s )) + 2A(R) 



T(s,s') = (d s T (s) ■ d 8 T (s)) + DkBT e V**T' . m 
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Finally we have to address the question of the reference state. Assuming with 
sufficient generality that the reference configuration of the chain is a straight line, 
thus Tq(s) = (se, where e is a constant vector [12] • At this point we can investigate 
also the saddle point of H[r (s)} with respect to (, the new variable, which reduces 
to the following two equations 

On general grounds ( can be non-zero only at zero temperature, when there are no 
fluctuations and the chain can indeed exhibit straight configurations. For any finite 
temperature straight configurations are not feasible without external constraints 
and we should rather have ( = 0. Taking this into account and solving the first of 
the above equations i.e. we obtain 

1 = (40) 
Let us now introduce a new parameter 



AK^ 

« - W < 41 > 

In D = A £ would be simply equal the renormalized persistence length. The factor 
4/3 obtained for D = 3 is a consequence of the 1/D ansatz, i.e. is a consequence 
of the approximate nature of the evaluation of the partition function. We do 
not attribute much importance to this discrepancy. With the introduction of £ 
equations Eqs. EE] can be reduced to a rather tame set of formulas 

B(a) = 2£( S + £(e-^-l)) 

T(s) = e~ s/i . (42) 

These two formulas deserve some interpretation. The first of the above two equa- 
tions obviously represents the average size of the chain of contour length s — s' 
since at the saddle point B{s — s') = <(r(s) — r(s')) 2 >. It is equal to the usual 
Kratky - Porod expression for the average end-end distance squared of a semiflex- 
ible chain. The second equation simply expresses the fact, that for a semiflexi- 
ble chain the orientational correlation function, which at the saddle point equals 
T(s — s') = <d s r(s) ■ d s tr(s')>, is exponentially decaying with a characteristic 
length equal to the persistence length. This again is completely consistent with 
the Kratky - Porod result for a semiflexible chain ^ 
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5 Renormalized bending rigidity 



The above developments now allow us to analyze the explicit form of the bending 
rigidity renormalization. Taking first of all into account the definition of B(s, s') 
and T(s, s') we can write the dipolar interaction potential on the saddle point 
diagonalized level in the following form 

V(B(s, s'), T(s, s')) = V R (B(s, s>))T(s, s') - V P (B(s, s'))B(s, s')T(s, s% (43) 

allowing us to evaluate 55Kc as a function of the parameters of this interaction 
potential. We obtain 

SK C = ~\jds s 2 (V R (B(s)) - V P (B(s))B(s)) - 

It is instructive to compare this with the analogous formulation of the problem for 
monopolar interactions [T2]. There it was derived that 

° 12 J dB{s) V ; 

In view of this, it seems that a possible interpretation of Eq. El would be that the 
orientational part of the interaction potential averaged over chain conformational 
fluctuations leads to effective intersegment attractions with a repulsive as well as 
attractive components. 

An exact evaluation of the highly non-linear Eq. EJ note that SKq is on 
the l.h.s. of this equation as well as on the r.h.s., hidden in £ of the definitions 
Eq. H21 is not feasible and we have to investigate the properties of its solution 
numerically. First of all we break the evaluation of the integral Eq. El into two 
parts SKc = SKq + 8Kq , defined as 

6Kg> = -- J dss 2 (V R {Ky/Bi^)-V P {Ky/B(sj)B( S : 
SKb = - T2 J a dss [-^^ V P {nVW)) ~ m{s) B{s 



where we have taken into account that the interaction potential can be written 
as a function of the argument «|r(s) — r(s')|. Furthermore the lower cutoff was 
set equal to a, of the order of the thickness of the chain, whose numerical value 
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was taken as a = lnm. This cutoff stemms from the breakdown of the continuum 
elasticity at small lengthscales. Contrary to the monopolar case ^2], this cutoff 
is essential and reflects the faster decay of the dipolar interactions with respect to 
the separation compared to the monopolar case. The next step is to write above 
relations in the form suitable for numerical evaluation. Let us first of all introduce 




Figure 1: The ( + F(() Eq. EZIfor four different values of na, viz. Ka = 100, 10, 1 
and 0.1 (from left to right). The dimensionless strength of the dipolar interaction 

(k£b) [ — ) was taken as I. At large values of £ the dependence is obviously linear. 



the following notation for the interaction potential 

2 KT 

V R (r) = /° 6 3 (l + Kr)=v K 3 f R ( K r) 
47ree r 6 

2 — KT 

V P (r) = -^^-(3 + 3Kr + (Kr) 2 ))=v K 5 fp(Kr), (47) 
47reeo r° 

where /r,p(w) are obviously defined as 

u 6 



f P (u) = —(3 + 3u + u 2 ). (4E 



u 



Introducing furthermore 



B(z) = z - 1 + e~ z with z = - (49) 
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and setting 

£=1/(>/2k£), (50) 
we can derive the following form for SK^ 

= -|«)3 £ dzz * U {J^rdj - f P (y^)/e) ^) ■ (si) 



and for 8Kq 



(52) 

where the derivative in and /p stands for /' = d ^^ > and thus we have 

e~ u 

f' R {u) = -— 5 (3 + 3u + u 2 ) 

e~ u 

/p(«) = -^( 15 + 15m + 6m 2 + m 3 ). (53) 

We should again note here that the equations for the renormalized bending mod- 
ulus Eqs. |^ and are essentially non-linear, since 8Kq is also hidden in the 
variable £ = dk B T / (A\/2k(K c + 5K C ))- Let us finally introduce a new variable ( 
defined as _ 

c= r = ^ = < 4 ^vi^, ,54) 

S S ^ 3fc B T 3 Ad V ; 

where Ad = 1/k is the Debye screening length. £ is nothing but the inverse reduced 
screeing length as introduced by Everaers et al. Furthermore if 

uKbT dksT 

then the definition of the renormalized bending rigidity K c = Kc + SKc can be 
cast into the form 

C = C (0) - (> (0 + g^ (2) (C)) • (56) 

Because of Eq. El we can furthermore write in 3 dimensions 

C + F(0 = c + J (f^) 2 C 3 (^ (1) (C) + £C 2jr(2) (0) = C (0) - (57) 
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Figure 2: Renormalized persistence length Eq. for three different values of 
= 20,10 and lnm (upper, middle, lower curves), circles. Also shown are the 
approximate solutions obtained from the expansion Eq. ESI The strength of the 

dipolar interaction £g (<^) was taken as 1 nm. 



This equation gives the functional dependence of ( on the parameters of the dipolar 
interaction, most notably \p. From the definition of ( Eq. EH the solution of the 
above equation immediately leads to the renormalized value of the persistence 
length. The functions F^(Q have obviously been defined as 

^ (2) (o = jT a ^/(^(cv^)-/p(cv^)-/p(cv / ^)c 2 ^) 



Fig. ^ shows the dependence of ( + F(() on £ for four different values of (na) that 
enters the definition of F(() via the lower bound of the integrals Eq. EH1 The 
solution of Eq. |37| is obtained from Fig. ^ by simply looking at the intersection 
between the curves on the figure with a line parallel to the ( axis at A more 
transparent form of these equations that will be solved numerically can be obtained 
if we go back to the original variables and can thus write for the renormalized 
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Figure 3: Fit (bold line)) of the renormalized persistence length to a power of 
X a for three different values of — 1, 10 and 20 nm (lower, middle and upper 
circles). The scaling exponent obtained from the fit is a = 0.75. The strength of 

the dipolar interaction l B ( — ) was taken as 1 nm. 



persistence length 



p p 4v/2 W ^3 A D J ^ ^3 / \ 3 A D J v 3 A D 

(59) 

This equation is the main result of our paper. Its solution gives the renormalized 
value of the bending rigidity or equivalently the persistence length as a function of 
the dipolar interaction parameters. In general it can only be solved numerically but 
simplified analytic solutions can be obtained in limited regions of the paramater 
phase space. 

It is instructive to compare these results with those derived for a polyelectrolyte 
chain with simple screened Debye - Hiickel monopolar electrostatic interactions 
along the chain, where the fixed charges are at separation A along the chain. If 
we formulate the results derived in in the same language as used above, we get 

< = C<0) + IFfe C2F(0 ' (60) 
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where in this case 



f°° f yC= — \ e c \ B(z) 
F(Q= dzz* 1 + C^B(z) _ V 3 . (61) 
■/o V / B^) 

The analysis of these equations (though written in a somewhat different, yet com- 
pletely equivalent form) has already been performed before and we direct the 
reader to that work 



6 Results 

We noted already that the basic equations derived for the renormalization of the 
persistence length of a semiflexible dipolar chain have no simple analytic solution. 
Figures El and 01 thus present numerical solutions of Eq. EHl for some values of 

parameters Ad, while we always consider the case y^j — ^nm~ l . This 
last parameter sets the overall scale of the renormalization and is not crucially 
significant so we do not study its variation in detail. Two universal behaviors 
seem to emerge. For small values of the screenenig length Ad the behavior is given 
in Fig. El For large values of the screening length an altogether different type of 
behavior is seen Fig. El Numerical solutions of Eq. EHldo not show any indications 
of a possible phase transition to an ordered phase at non-zero temperatures though 
due to the attractive component of the dipolar interactions one might expect them 
to. This is completely consistent with the behavior of the flexible chain jH]. 

Let us investigate if we can understand the two types of behavior depicted on 
Figs. El and El with simple analytical arguments. What we basically need to figure 
out is what would be the typical contributions to the integrals Eq. EElthat enter 
the equation for the renormalized bending rigidity Eq. 023 It is instructive to 
consider first the case of small screening lengths. Here the main contribution to 
the integrals Eq. EH1 along the polyelectrolyte comes from short length scales. In 
this limit one derives that B(s) ~ s 2 and locally the chain thus behaves in the 
Kratky - Porod style ^0j- This leads to the following approximate form of F(() 

F(() _ , (, W(C) + g F(2)(c) ) . 2 ^J2 )k f (-» - 18 i| - . (^|) 2 + (i|) 3 ) . 

(62) 

Eq. 0211 in this very same limit then assumes the form 

L „ = Lf ^ (-25 -7 (m) + («,)>-'"» _ 3 Bi((M) ) (63) 
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where E\(u) is the standard exponential integral with a limiting form of lim u >o E\{u) = 

—7 — log u + 0(u). Except for very small values of the screening length where obvi- 
ously L P R ^ m L p °\ a valid approximation of Eq. EB1 capturing its dominant features 
would be 

^4°' + ^fe) 2 log^. (64, 

This result could be seen as the dipolar chain equivalent of the Odijk - Skolnick 
- Fixman [2j result for the monopolar chain. In that case the behavior of the 
persistence length at small values of the screening length is L p ~ X 2 D . The 
difference between this result and Eq. EH is purely due to the swifter decay of 
the dipolar interaction potential vs. separation, r, in monopolar case r" 1 and 
in dipolar case r -3 . This can be seen crudely as follows: since the interaction 
potential for the dipolar chain falls off twice as fast as the monopolar potential the 
OSF reasoning should give for the electrostatic persistence length X 2 D x \ D 2 = \° D . 
The zero in the exponent translates as usual into the log, which is exactly Eq. 

The way the approximate result Eq. EH] fares when compared with the full 
numerical solution of Eq. EH1 is shown on Fig. [21 Obviously the approximaton 
is not bad for small enough screening length but becomes progressively worse for 
larger bare persistence lengths and larger screening. The same reasoning applied to 
the general multipole of order n would lead to the electrostatic persistence length 
with scaling A^ n . One can thus hardly expect any electrostatic renormalisation 
effects above dipolar. 

Next we consider the case of large(er) screening lengths. Here the main con- 
tribution to the integrals Eq. EB1 along the polyelectrolyte comes from large(er) 
length scales. In this case the local properties of the chain are characterized by an 
essentially free flight behavior 10 j of the chain, leading to B(s) ~ 2£s. Here we 
can derive 

F(C) = C 3 (W) + ^ (2) (C)) = - dzz" 2 (2 + 2C^z + Cz) = ~C 3 , 

(65) 

which leads to the following approximate form of Eq. EH 

(66) 

This equation can have two possible solutions depending on the magnitude of L p °\ 
First of all if Lp 3> L p °\ meaning that the value of the persistence length is 
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determined mostly by the electrostatic interactions along the chain, we have 




2 



) 



1/4 




(67) 



This approximate result is nicely exhibited also on Fig. El where it is seen how the 
exact numerical solutions of Eq. EH] approach the above scaling limit for sufficiently 
large values of the screening length with the scaling exponent exactly equal to the 
above prediction. The form Eq. EH is thus a limiting law for the electrostatic part 
of the persistence length valid universally for large screening lengths. One should 
be however aware, that in the above calculation by assumption the length of the 
chain is always the largest length in the system. 

In the other case with a predominant contribution of the bare persistence 
length, i.e. L V R ^ ~ thus in the limit where the strength of the dipolar in- 
teraction is very small, we end up with the following limiting form 



The latter is only valid for small enough strength of the dipolar interaction and is 
difficult to discern in the numerical solution. Since all our numerical investigations 
on Figs. 01 are based on the assumption of a fairly large strength of the dipolar 
interactions, this regime is not distinguishable on the graphs. 

These last two results could again be compared with the expressions valid for 
the monopolar chain in the same regions of the parameter space and derived within 
the same theoretical framework [T2]. In that case the behavior of the persistence 
length at large values of the screening length would be L P R ^ ~ A^>, where (3 is 
either 7/6 or 7, depending on whether Lp / L p ^ is either large or small. Again the 
difference between these results and Eqs. EH and EB1 is purely due to the swifter 
decay of the dipolar interaction potential vs. separation; monopolar case r -1 and 



into the second term of Eq. EE1 wherefrom the exponents (3 = 3/4 or 3 can be 
derived straightforwardly. It is thus clear that the behavior of the electrostatic 
part of the persistence length of a monopolar and dipolar chain are intimately 
related in both relevant limits of weak and strong screening. Applying again the 
same reasoning to the general multipole of order n would lead to the electrostatic 
persistence length with scaling (3 = (7 — 2n) or (3 = (7 — 2n)/(6 — n). Again 
electrostatic renormalisation effects above dipolar are very small. 

The remaining question in this context would be how robust are the two regimes 
derived above that lead to the approximate forms Eqs. EH and 16 71 '' In assessing 




3 



(68) 



dipolar case r 



3 . This difference introduces an additional factor of A 
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the range of validity of these different approximations we can invoke a related sit- 
uation already encountered in the case of a monopolar chain ^2]- In that system 
extensive simulations |S1 IH E] left no doubt that the OSF regime, corresponding 
in the context of the dipolar chain to Eq. |^ is very robust and extends over a 
broad region of the parameter phase space. The sub-OSF laws for the electrostatic 
renormalization of the persistence length giving f3 in the vicinity of 1 were effec- 
tively ruled out. Translating this into the context of the dipolar chain would make 
the range of validity fo Eq. EHlfahly narrow. But these are all conjectures since we 
are aware of no detailed simulations for semiflexible dipolar chains, setting aside 
the extensive work by Muthukumar on the flexible chain jH] , that would match the 
superb work perfomed recently in the context of the monopolar chain |H1SE]- 

7 Conclusions 

We presented an analysis of the electrostatic contribution to the persistence length 
of a semiflexible screened-dipolar chain. The formal context of our analysis is pro- 
vided by the l/Z)-expansion method that has been already successfully applied to 
different problems of polymer physics [121 HH1 HHj • 1 /-D-expansion is closely related 
to different variational schemes that have been amply applied to the problem of 
electrostatic rigidity of charged polymers (201 |2U 1221 12H 125- What singles out 
our approach is the fact that we work consistently with a semiflexible Hamiltonian 
and enforce the inextensibility constraint on a global level. This leads to some dis- 
crepancies between variational formulations and 1 / D-expansion method. However 
the robustness of the OSF regime transpires quite clearly from the 1/D-expansion 
in the context of a monopolar polyelectrolyte, giving us some confidence that an 
analogous result for a dipolar chain would have the same range of validity. 

The main step in our formulism was to find an appropriate way to treat the 
orientational dependence of the intersegment interaction potential. This has been 
accomplished by introducing additional auxiliary fields that in their turn lead to 
new saddle-point equations. Though the derivation of the renormalized elasticity 
is as a consequence more complicated, it leads to manageable and transparent 
results. The main conclusion based on these results would be that the dependence 
of the electrostatic persistence length on the screening length for a dipolar chain 
is much less pronounced then for a monopolar chain. This is clearly seen in the 
case of large as well as small screening. Also it appears that a semiflexible chain 
does not give rise to any localized structures as those described by Muthukumar 
[Hj in the context of a flexible chain. This is indeed not surprising since the chain 
elasticity is obviously strong enough to prevent extensive looping of the chain that 
would lead to local aggregation of the chain based on the attractive component of 
the dipolar interaction. This attractive component is also not strong enough to 
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lead to any type of phase transition to an ordered state at non-zero temperatures. 
This conclusion is just as valid for a semiflexible as it is for a flexible chain. 

Eventually the conclusions arrived at in this work would have to be tested 
against extensive simulations just as they were in the case of a monopolar charged 
chain. As for the experimental situation where repeated claims of a sub-OSF 
regime have been voiced [23 1213, one wonders if they could in fact be a conse- 
quence of a specific association of the mobile counterions with fixed charges on the 
polyelectrolyte backbone producing a system not far from the model dipolar chain 
analyzed in this work. In order to test this hypothesis one would have to estimate 
independently the effective dipole moments of the chain segments. The results 
presented in this work could then serve as a guideline to differentiate between 
monopolar and dipolar OSF-like behavior. 
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